Search for millicharged particles in reactor neutrino experiments: 
a probe of the PVLAS anomaly 
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It has been recently suggested that the vacuum magnetic dichroism observed by the PVLAS 
experiment could be explained by pair production of a new light, m e ~ 0.1 eV, millicharged, 
q e ~ 3 • 10 _6 e, fermions (e). In addition, it has been pointed out that millicharged particles with 
q e > 10 _9 e appear naturally in models based on the string theory. 

We show that low energy reactor neutrino experiments provide a sensitive probe of millicharged 
particles. Considering, as an example, recent results of the TEXONO experiment searching for 
neutrino magnetic moment, a new upper bound q E < 10~ 5 e for the mass region m s < 1 keV is 
derived. These results enhance motivations for a more sensitive search for such particles in near 
future experiments. Furthemore, a direct experimental limit on the electric charge of the electron 
antineutrino g Sc < 3.7 ■ 10 _12 e is obtained. 

PACS numbers: 14.80.-j, 12.20.Fv, 13.20.Cz 
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Several realistic extensions of the standard model 
of elementary particle physics suggest the existence of 
particles with small, unquantized electric charge. They 
appear naturally and generically from kinetic- mixing in 
theories containing at least two U(l) gauge factors, e.g. 
in paraphoton models By adding a second, hidden 
photon (the "shadow" photon) one could construct grand 
unified models which contain particles with an electric 
charge very small compared to the electron charge Q. 
These considerations have stimulated new theoretical 
works and experimental tests, see e.g. [1, 0, H, 0] and 
references therein. 

Recently, a laboratory experiment PVLAS has re- 
ported on observation of an anomalous rotation of the 
polarization plane of a photon beam incident on a region 
of a static transverse magnetic field in vacuum Q . The 
PVLAS signal appears to be in contradiction with the 
negative results of the CAST experiment on search for 
solar axions [§] and astrophysical constrains [l(| (how- 
ever, one may argue that astrophysical constraints are 
also model dependent and have considerable theoretical 
and experimental uncertainties). 

Although many models explaining the PVLAS-CAST 
discrepancy appeared recently , not all of them allow 
for experimental tests. One of the experimentally inter- 
esting explanation suggests that the PVLAS anomalously 
large effect may originate from pair production of light, 
m £ ~ 0.1 eV, millicharged, q £ ~ 3 ■ 10 _6 e, fermions (e), 
whose existence is not excluded by the present experi- 
mental data [ll|. It is also pointed out that the multi- 
ple U(l) factors, the size of kinetic-mixing, and suitable 
matter representations required to explain the PVLAS 
data occur very naturally in the context of realistic ex- 
tensions of the Standard Model based on string theory 
[l3|. The string scale ~ 10 11 GeV predicts the existence 
of particles with q e > lCP 9 e, thus making the region 
g £ ~ 10~ 9 — 10 _5 e of great interest for searching for mil- 



licharged particles in future experiments. 

In this note, we show that reactor neutrino experi- 
ments searching for non-standard neutrino interactions 
provide a sensitive probe to light millicharged particles. 
We discuss possibilities for more sensitive searches for 
such particles in near future experiments. 

It is well known that since the electromagnetic cross 
section is orders of magnitude larger than the weak cross 
section, the existence of non-zero magnetic moment of 
the electron antineutrino would increase the total rate of 
events in v — e — - scattering experiments [hjj . The in- 
tegral cross section for neutrino magnetic scattering 
depends very weakly on the neutrino energy. It rises only 
logarithmically with the neutrino energy, while the total 
neutrino cross section rises linearly with neutrino energy 
E v . Thus, it is advantageous to search for neutrino mag- 
netic moment using low neutrino energy. These consid- 
erations point to nuclear reactors as probably the most 
convenient source of neutrinos. Reactors with power of 
2 GW emit about 4 x 10 20 v e per second [l5j]. The emis- 
sion energy spectrum is quite well understood. It has a 
broad distribution over energies up to 10 MeV, with a 
maximum intensity at 0.5-1.5 MeV [bj ]. 

Consider, as an example, recent results of the TEX- 
ONO collaboration searching for electron anti-neutrino 
magnetic moment 16]. The experiment has been set 
up at the Kuo-Sheng Nuclear Power Station at a dis- 
tance of 28 m from the 2.9 GW reactor core. A detector 
threshold of 5 keV and a background level of 1 counts 
day _1 keV _1 kg _1 at 12-60 keV was achieved with a high 
purity germanium detector of mass 1.06 kg surrounded by 
anti-Compton detectors with Nal(Tl) and CsI(Tl) crystal 
scintillators. 

The excess AiV„_ e of events in the Ge detector due to 
the v — e magnetic scattering is given by 
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where E e is the electron recoil energy, dN v jdE v is the 
reactor neutrino flux, da^jdE,, is the v — e scattering 
differential cross section via the v magnetic moment, and 
k is a factor corresponding to the convolution of the Ge 
detector acceptance and mass. 

The cross section da^/dE e is given by, see e.g. 
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where p v is the neutrino magnetic moment and /is = 
e/2m e is Bohr magneton. 

One can see that for the electron energy interval from 
E rnin = 12 keV to E™ ax = 60 keV, analyzed by TEX- 
ONO [HI, the cross section of v e — e scattering due to 
non-zero magnetic moment is practically independent of 
energy for E v > 0.5 MeV and is equal: 



a u = 4.3 x 10 



-45 



fJ'U 



10" 



j Pb 



cm 



(3) 



Thus, the total number AN v - e of events collected during 
period of time At could be estimated as 



AN^ e = Ua^NeAt 



(4) 



where /„ = 5.8 x 10 12 cm" 2 s" 1 is the integrated reactor 
neutrino flux ( the detection efficiency of the Ge detector 
is ~ 100% for energy deposition above 10 keV) [l(| and 
N e is the total number of electrons in the detector. 

Using the above values we found that the TEX- 
ONO experiment should detect an excess of AiV„_ e ~ 
5.8 events/day if [i Ve = 2 x 10 _1 Vb- This value is m an 
excellent agreement with the number of events quoted in 
ref. [3. 

Having this in mind, consider the case of millicharged 
particles. The nuclear reactors could also be a source of 
such light particles. Indeed, the 2.9 GW power reactor 
core emits large number of 7-quanta, more than 5 x 10 20 
per second, and these gammas could convert within the 
reactor into pairs of light e-particles. Even a small frac- 
tion of these particles could lead to an observable excess 
of electrons from e — e" scattering in the detector whose 
experimental signature is identical to that of v — e~ mag- 
netic scattering events. 

Consider now the number of e — e events in the TEX- 
ONO Ge detector. The excess AiV e _ e of events due to 
e — e scattering in the detector is given by: 
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where dN £ /dE e is the assumed e-particles flux, and 
da E /dE e is the e — e scattering differential cross section. 

For E E ^> E" lm ^> m e , i.e. m e < 1 keV, the cross 
section da e / dE e is independent of E e and is given by 
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where t-q is the classical electron radius, or 



a F = 2x 10" 



(7) 



= 12 keV to 



for the recoil energy region from E r e 

E rnax = qq ke y 

The highest rate mechanism for the light e-particle pro- 
duction by ~1 MeV gammas is the ee-pair production in 
the Compton scattering 
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The cross section (a ££ ) of this process is proportional 
to q £ and is characterized by q £ , m £ and spin of the e- 
particles. It dominates over Bethe-Heitler pair produc- 
tion which is proportional to q £ . 

For the given photon energy interval (E 7 ; E 1 + dE 7 ) 
the total number of e-particles per second produced in 
the reactor core can be calculated as: 
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where dN-y/dE-y is the reactor gamma spectrum and L is 
the mean free path of photons in the reactor core given 
by 
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where <r tot , N A , A, Z, p are the total photon cross section, 
the Avogadro's number, the atomic weight and charge, 
and the density of the core material (we assume for the 
moment that this material has uniform composition). 
The cross section a F£ is found to be 
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where ac is the Compton scattering cross section. It is 
calculated under assumption a ££ -C <7c an d is valid for 
a/3n(q £ /e) 2 ln(E 1 /2m £ ) < 1. 

Although the reactor gamma spectrum has all possible 
energies up to 10 MeV, not all photons effectively con- 
tribute to the ee pair production through the reaction of 
Eq.flSJ). Indeed, for high Z-materials, Z ~ 90 - 100, in 
the low energy region E-y < 1 MeV the total gamma cross 
section is dominated by the atomic photoelectric effect, 
while at energies higher than E 1 > 5 MeV the e + e~ pair 
production dominates. Thus , the suitable gamma energy 
interval for the effective ee pairs production via reaction 
of Eq. © is 1 - 5 MeV. For gammas with these ener- 
gies the total cross section is dominated by the Compton 
scattering for any material with Z up to 100. Hence, 
we may assume a tot — &c- Combining this with Eq.©, 
Eq. (fT0"l) and Eq. lfTTj) one can see that the total number 
of e's produced inside the reactor per sec is 
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Thus, in this approximation the e-flux is independent on 
fuel composition, reactor core geometry, etc. 

The total number of photons produced in reactor core 
can be estimated by using results on 7-energy spectra 
measurements [l7j . used in the experiment searchin g fo r 
axion decays at a nuclear power reactor at Bugey [18j |. 
The reactor gamma spectrum receives contributions from 
several processes following the nuclear fission process, see 
e.g. [HI]. The comparison of numbers of neutrinos [lij 
and gammas [l?} emitted per fission for various energy 
intervals shows that the total number of gammas emit- 
ted in the energy range 1-5 MeV for the nuclear fuels 
of interest, i.e. 235 U, 239 Pu, 238 U, 241 Pu, and 252 Cf, is 
comparable with the corresponding number of neutrinos. 
For example, for 235 U the number of neutrinos and gam- 
mas emitted in the 1-5 MeV energy range is found to be 
N u = 5.7 and N 1 — 2.1 per fission, respectively. The ra- 
tio of the number of neutrinos to the number of gammas 
are the same within ±15% for different type of fuels. For 
the energy range 0.15 - 0.5 MeV the number of gammas 
is about 30-40% higher than the number of neutrinos. 
The energy interval 5-10 MeV is not very important, be- 
cause less than about 5% of particles have those energies. 
Conservatively, we may assume that N 1 = 0.3NI of the 
neutrino flux N' v in the energy interval 1- 5 MeV, or 

N 1 = 0.2N U (13) 

of the total neutrino flux N v . This is calculated by us- 
ing reactor neutrino spectra for different nuclear fuels 
described in ref. 

Similar to Eq.([4]) the total number AiV e _ e of e — e 
events collected during period of time At is 



AiV £ _ e = f e a e N e At 



(14) 



where f £ is the flux of e's with energy E e ^> £™ m . The 
energy of e's produced by 1- 5 MeV gammas is found to 
be widely distributed over the range < E e < 5 MeV 
with a few %'s of events lying in the low energy region 
E £ < 50 keV. 

Combining Eq.{l2]) and Eq. (fT3| . the flux of e's per sec 

is 
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Finally, using Eq.© and Eq. ([T5|) and taking /„ = 5.8 x 
10 12 cm" 2 s- 1 and N e = 2.8 x 10 26 0, one can find 
that the number of e — e events expected to be seen in 
the TEXONO Ge detector per day is 
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Comparing this to the upper bound on the number of 
events implied by the upper bound on the neutrino mag- 
netic moment fx v < 1.3T0~ 10 /xb [l6| and assuming E^ — 
1 MeV we obtain : 



for the mass region m £ <1 keV. This estimate is weakly 
depended on the e-mass, variations of the selected recoil 
energy interval and uncertainties of the e-flux. It gives 
an order of magnitude for the upper limit on q e and may 
be strengthened and extended to the mass region m e > 
1 keV by more accurate and detailed Monte Carlo sim- 
ulations. The limit of Eq. (fT7|) is valid for a e-particle 
lifetime r e [s] > 10 _13 m e [eV]. The attenuation of the flux 
due to (electromagnetic) interactions in the shielding is 
found to be negligible, since for the limit of Eq. (fl7| the 
e-particle mean free path in iron is ^>1 km, as compared 
with the total shield thickness of < 28 m used at the 
Kuo-Sheng Nuclear Power Station. 

The results obtained could be also used to derive direct 
experimental limit on the electron anti-neutrino electric 
charge: 



< 3.7 x 10 
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which is within two orders of magnitude of astrophysical 
bounds 

The statistical limit on the sensitivity of a reactor ex- 
periment searching for millicharged particles is propor- 
tional to g 4 and is set by the value of q e . Thus, to im- 
prove the limit of Eq. (|17p large target mass is required. 
Because of the 1/E e dependence of the e — e scattering 
cross section the detection of low recoil electron energies 
is also crucial for the improvement of sensitivity of further 
searches. A possible approach, which could satisfy both 
requirements, is to use a massive liquid argon detector. 
An example of such detector is the one developing in the 
framework of the ArDM project for the search of Dark 
Matter (20| . This experiment aims at running almost 
1 ton of effective mass detector and is capable to mea- 
sure deposited energies down to 10 keV. The significance 
of the e-particles discovery with such detector, assuming 
q e = 3 • 10 _6 e, scales as 



s*i<r>(™) 1/a 



(19) 



q e < 10" 5 e 
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where S = 2(-s/n~ s + rib ~ \fn~b) [2l[ and n s = 10~ 2 Mt, 
nt = noMt are the number of detected signal and back- 
ground events, respectively. Here M is the mass of the 
detector in kg, t is the running time in days and no is the 
background counting rate in kg _1 day _1 for the energy in- 
terval of interest. Thus, assuming S > 3, M ~ 10 3 kg and 
t ~ 10 2 days, a background level of no < 1 event/kg/day 
has to be achieved [2(| ■ The experiment could also signif- 
icantly improve the recently obtained bounds on neutrino 
magnetic moment [l6j |. [22j |. [23| pushing it down to the 
region [i v ~ 10~ 12 /ib. A bound in this region is com- 
parable with the astrophysical limits [19j and will be of 
great interest for many extensions of the Standard Model 
[24| . Note, that e — e events could be distinguished from 
the v — e magnetic scattering due to the (1/E e ) 2 profile 
of their recoil energy spectrum, while v — e events have 
1/E e energy distribution. 

Another experimental approach relies on the precise 
study of the positronium atom 25]. If low-mass mil- 
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licharged particles exist, the orthopositronium (o — Ps), 
the triplet bound states of electron and positron, could 
decay " invisibly" into ee pairs since the e-particles would 
penetrate any type of calorimeter mostly without inter- 
action [26|, [23] . For spin 4 millicharged £-particle and for 
m e <C m e the best experimental bound obtained by the 
ETH-INR experiment searching for o — Ps — » invisible 
is q e < 3.4 • 10 -5 [28|, which is comparable with the limit 
ofEq.CTJ. 

The sensitivity of experiments on invisible decay of or- 
thopositronium is proportional to , thus making further 
searches of this decay mode more interesting and promis- 
ing [25j]. The ETH-INR experiment achieved sensiti vity 
to Br(o-Ps -> ee) as small as 4.2 • 1CT 7 (90% C.L.)[lf. 
To intersect the PVLAS parameter space a new experi- 
ment with the sensitivity improved by a factor of 100 is 
required. 



Recently, a method based on the measuring of the elec- 
tric current comprised of the millicharged particles pro- 
duced in a cavity has been proposed [2!|, see also [30]. 
The millicharged particles, if they exist, contribute to 
the energy loss of the cavity and thus leave an imprint 
on the cavity's quality factor. Already conservative esti- 
mates performed in [29j show that this experiment will 
substantially bound the parameter space for e-particles 
allowing firmly exclude ( or confirm) the region suggested 
by the PVLAS signal. 
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